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ABSTRACT 

Using CGS4 on UKIRT, we have obtained 2.00-2.45 /im infrared spectra of a number 
of low-mass X-ray binaries including Sco X-l, Sco X-2, and GX13+1. Sco X-l shows 
emission lines only, supporting our previous conclusion that the spectral type of the 
evolved secondary must be earlier than G5. Emission lines are also seen in the spectrum 
of Sco X-2, confirming the identity of the IR counterpart. We report the detection of 
CO bands in GX13+1 and estimate the most likely spectral type of the secondary 
to be K5iii. We also find P Cygni type profiles in the Brackett 7 lines of Sco X-l 
and GX13+1, indicating the presence of high velocity outflows in these systems. We 
present spectra of candidate IR counterparts to several other elusive X-ray binaries. 
Finally, implications for the nature and classification of these systems are discussed. 

Key words: infrared: stars - X-rays: stars - binaries: close - binaries: spectroscopic 
- accretion, accretion discs 



1 INTRODUCTION 

In low-mass X-ray binaries (LMXBs), a neutron star or 
black hole primary accretes material from its late-type com- 
panion star, producing intense X-ray emission. LMXBs can 
be divided into subclasses according to location within the 
Galaxy, accretion characteristics, and luminosity (e.g. van 

' Paradijs & McClintock 1995). The "bright Galactic bulge" 
sources (GBS) are located within 15° longitude and 2° lat- 
itude of the Galactic Centre (see e.g. Warwick et aJ. 1988) 

' and are among the most luminous X-ray sources in the 
Galaxy (typical L x ~ 10 38 ergs -1 ). The GBS have shown 
no X-ray bursts and attempts to detect orbital variability 
have been unsuccessful, suggesting that their periods may 
be longer than those of canonical LMXBs (Charles & Naylor 
1992, hereafter CN92). In addition, heavy obscuration in the 
Galactic Centre region has made optical study of the GBS 
nearly impossible. There are as yet no observations which 
explain the dichotomy between the poorly understood GBS 
and the rest of the LMXBs. The most likely theories suggest 
that the secondary stars in the GBS are late-type giants, in 
contrast to the quasi-main sequence stars in other LMXBs. 

However, the infrared provides us with an ideal window 
for observing these systems. The IR has two primary ad- 
vantages: the late-type secondaries in LMXBs are brighter 
relative to the accretion discs, and, more importantly for 
the GBS, the ratio of V- to if-band extinction is nearly 
10 (Naylor, Charles, & Longmore 1991, hereafter NCL91). 
Over the past eight years, we have developed a program of 
IR observations of X-ray binaries (XRBs), beginning with 



the discovery via colours or variability of candidates for the 
IR counterparts to the X-ray sources using precise X-ray and 
radio locations (NCL91). Following this photometric survey 
we began an IR spectroscopic survey of LMXBs. In 1995 
we obtained IR spectra of the LMXB Sco X-l and the GBS 
systems GX1+4 and GX13+1; these results were published 
in Bandyopadhyay et aJ. (1997), hereafter Paper I. 

Continuing our spectroscopic survey of LMXBs, in this 
paper we present new results obtained from UKIRT in 1997. 
In Paper I, we presented a JHK spectrum of the prototype 
LMXB Sco X-l taken in 1992 with an older, low resolution 
array and a short integration time. We have now observed 
Sco X-l at higher resolution and with a substantially longer 
exposure time. We have also obtained the first Jf-band spec- 
trum of the IR counterpart to the GBS Sco X-2 (GX 349+2), 
which was initially identified with a variable radio source 
(Cooke & Ponman 1991); variability in the _R-band coun- 
terpart was later found by Wachter & Margon (1996). In 
Paper I we presented a spectrum of the heavily obscured 
GBS GX13+1 which confirmed the identity of the IR coun- 
terpart. We now present a new spectrum of this source which 
clearly shows the CO bands and metal lines of its late-type 
secondary. Finally, we have obtained IR spectra of candidate 
stars within the X-ray/radio error circles of the GBS GX5-1, 
GX17+2, and the LMXB 4U2129+47. 
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Table 1. Journal of observations 



Object 


Date 


UT 

(hrs) 


Exposure time 
(sees) 


Sco X-l 


1/7/1997 


6:55 


3200 


Sco X-2 


1/7/1997 


8:35 


2880 




1/7/1997 


9:38 


1920 


GX13+1 


1/7/1997 


10:48 


2880 




1/7/1997 


11:55 


2400 


GX17+2 


2/7/1997 


8:50 


2400 


GX5-l a #502 


2/7/1997 


10:10 


960 


GX5-l a #503 


2/7/1997 


10:49 


960 


4U2129+47 


3/7/1997 


12:35 


3360 



a Stars 502 and 503 are candidate stars within the X-ray error 
circle for GX5-1. 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained A'-band (2.00-2.45 fim) spectra using the 
Cooled Grating Spectrometer (CGS4) on the 3.8-m United 
Kingdom Infrared Telescope on Mauna Kea during the 
nights of 1997 July 1-3 UT. The 75 1/mm grating was used 
with the 150 mm camera and the 256x256 pixel InSb ar- 
ray. Target observations were bracketed by observations of 
A-type stars for removal of telluric atmospheric features. A 
journal of observations is presented in Table 1. 

The standard procedure of oversampling was used to 
minimise the effects of bad pixels (Wright 1995). The spectra 
were sampled over two pixels by mechanically shifting the 
array in 0.5 pixel steps in the dispersion direction, giving a 
full width half maximum resolution of 34 A (~460 kms _1 at 
2.25 lira). We employed the non-destructive readout mode 
of the detector in order to reduce the readout noise. The 
slit width was 1.23 arcseconds which corresponds to 1 pixel 
on the detector. In order to compensate for the fluctuating 
atmospheric emission lines we took relatively short expo- 
sures and nodded the telescope so that the object spectrum 
switched between two different spatial positions on the de- 
tector. For some objects, the slit was rotated from its default 
north-south position to avoid contamination of the target 
spectrum by nearby stars. Details of the design and use of 
CGS4 can be found in Mountain et al. (1990). 

The CGS4 data reduction system performs the initial 
reduction of the 2D images. These steps include the appli- 
cation of the bad pixel mask, bias and dark subtraction, 
flat field division, interlacing integrations taken at different 
detector positions, and co-adding and subtracting the nod- 
ded images (see Daly & Beard 1994). Extraction of the ID 
spectra, wavelength calibration, and removal of the telluric 
atmospheric features was then performed using IRAF. A 
more detailed description of the data reduction procedure is 
provided in Paper I. 



3 RESULTS 
3.1 Sco X-l 

Our K-band spectrum of Sco X-l is shown in Figure 1. 
Strong emission lines of HI, Hel, and Hell are the domi- 
nant features; no absorption lines from the secondary star 



Table 2. Measured line wavelengths and equivalent widths 



Line 


A 


Sco X-l 


Sco X-2 


GX13+1 




((jm) a 


(A) 


(A) 


(A) 


Hel 


2.059 


-2.8±0.4 


-26.0±1.6 


-3.7±0.6 


Hell 


2.098 


-1.0±0.2 






Hel 


2.114 


-1.0±0.2 






Br 7 


2.165 


-12.6±0.5 


-29. Oil. 8 


-4.8±0.4 


Hell 


2.187 


-2.7±0.4 






Cal (triplet) 


2.262 






1.2±0.5 


Mgl 


2.283 






0.8±0.4 


12 CO (2,0) 


2.295 






5.8±0.6 


12 CO (3,1) 


2.324 






3.2±0.6 


13 CO (2,0) 


2.345 






0.6±0.3 


12 CO (4,2) 


2.353 






3.7±0.6 


13 CO (3,1) 


2.377 






1.6±0.5 


12 CO (5,3) 


2.384 






3.8±0.7 


12 CO (6,4) 


2.416 






1.6±0.7 



a Typical wavelength calibration errors are approximately ±0.001 
(im. 



can be distinguished. A list of identified lines is found in 
Table 2. 

3.1.1 The secondary star 

In Paper I we presented a JHK spectrum of Sco X-l, taken 
with the earlier 64x64 CGS4 array. The spectrum showed no 
absorption features, and we concluded that the secondary is 
most likely a subgiant (see Paper I for a detailed discussion) . 
By modelling the appearance of two spectral template stars 
in the if- band at the distance and reddening of Sco X-l, we 
determined that the mass-donating star in Sco X-l is of a 
spectral type which shows little or no CO features, i.e. earlier 
than G5 for a subgiant (Kleinmann & Hall 1986, hereafter 
KH86). However, because the spectrum was taken with a 
short integration time on the older, lower resolution CGS4 
array, there remained some doubt about this conclusion. 

The Sco X-l spectrum presented here removes this 
doubt. With the new array and a ~53 minute integration 
time, there is still no evidence of absorption features in- 
trinsic to the secondary star in the K-b&nd spectrum, only 
the emission lines expected from the disc and/or the heated 
face of the secondary. The relatively long orbital period (P 
= 18. 9h) and the high mass transfer rate favour the for- 
mation of a large accretion disc (see Beall et al. 1984) and 
the optical spectrum is dominated by the disc (Schachter, 
Filippenko & Kahn 1989); consequently we might expect to 
see emission from the disc in the IR. However, based on our 
earlier modelling (see Paper I) , in our new spectrum it is un- 
likely that contamination by the disc would be sufficiently 
high to completely obscure the strong CO features expected 
in the late-type subgiants. Therefore the mass-donating star 
in Sco X-l is of at least an early G-type, and perhaps even 
earlier, making Sco X-l somewhat unusual. The majority of 
identified secondaries in LMXB systems have been of a later 
type, most frequently K or M stars. As such, it is interesting 
to note that Sco X-l, the "prototypical LMXB", may not 
be so typical after all. 
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Table 3. P Cygni line profile parameters a 





Sco X-l 


GX13+1 


Adopted rest velocity of star 


-230±40 


-70±30 


Red edge of Br 7 emission 


-420±30 


-430±30 


Blue edge of Br 7 absorption 


-2830±50 


-2480±40 


Wind outflow velocity 


-2600±60 


-2400±50 



All values given in km s 1 . Errors quoted are formal errors only. 



3.1.2 The P Cygni profile 

Despite the lack of absorption features from the secondary 
star, one absorption feature is present in our Sco X-l spec- 
trum. The Br 7 emission line has a noticeable absorption 
dip on its blue edge: a P Cygni profile, indicative of a wind 
in the system. By subtracting the measured velocity of the 
star from the velocity of the blue edge of the absorption fea- 
ture, the maximum outflow velocity of the wind is obtained. 
Although our spectrum has insufficient velocity resolution 
for an accurate measurement, we have made estimates of 
the radial velocity of Sco X-l (using the Hel and Hell emis- 
sion lines) and for the blue and red edges of the P Cygni 
profile; the results and the calculated outflow velocity are 
listed in Table 3. A P Cygni profile has also been observed 
in the LMXB GX1+4; the ~250kms _1 outflow in that sys- 
tem likely comes from the M5iii donor star (Chakrabarty, 
van Kerkwijk, & Larkin 1998). However, the wind velocity 
we infer for Sco X-l, approximately 2600 kms -1 , is super- 
sonic, and thus unlikely to originate from the companion 
star (e.g. Warner 1995, Dupree 1986). As in CVs, the wind 
probably originates from the hot accretion disc, which pro- 
duces a radiatively driven outflow. Given the dissimilar sizes 
of the two systems (P or b ^ 260d for GX1+4; Chakrabarty 
& Roche 1997) and the contrasting mechanisms for driving 
the wind, the order-of-magnitude velocity difference between 
the winds in GX1+4 and Sco X-l is not unexpected. Indeed, 
outflow velocities similar to that seen in Sco X-l have been 
obtained in models of disc winds in CVs (e.g. Mason et al. 
1995). However, to make accurate measurements of the out- 
flow velocity and obtain information about the origin of the 
wind in Sco X-l (e.g. from the shape of the line profile), high 
resolution observations of the Br 7 line are required. 

3.1.3 Shape of the continuum 

We have measured the slope of the Jf-band continuum of 
Sco X-l, after dereddening the spectrum using Eb-v = 0.3 
(determined from the 2200A interstellar absorption feature; 
Vrtilek et al. 1991). We fit the data with a power law of 
the form F\ oc \ a to represent the emission from a stan- 
dard steady state accretion disc (e.g. Frank, King, & Raine 
1992). Using this model, we find a power law index q= - 
3.17±0.02 for the IR spectrum of Sco X-l. Shahbaz et al. 
(1996) found that the optical spectrum of Sco X-l showed 
no evidence for any contribution from the secondary star 
and was well fitted with a— -2.46; however, they used an 
older estimate of the reddening (Eb-v = 0.15). Using a 
reddening of 0.3 and normalizing the flux to the observed 
V magnitude of Sco X-l (van Paradijs 1995) by folding the 
spectrum through the filter response, we find a — -2.78 for 



Table 4. Continuum slopes: power law indices (a) 





Sco X-l 


Steady State Disc 




Observed 


T out = 4200 K 


Optical 


-2.78±0.09 


-2.6±0.1 


IR 


-3.17±0.06 


-3.5±0.1 



the optical spectrum. We then estimated the expected outer 
disc temperature (T out ) for a steady state disc with no ir- 
radiation and found T out ~4200 K. The continuum slopes 
of the observed optical and IR spectra are different from 
that expected from a steady state disc at this temperature; 
these results are summarized in Table 4. Figure 2 shows the 
observed optical and K-b&nd spectra plotted with the sim- 
ulated steady state disc for comparison. 

What is the cause of these discrepancies? As Sco X-l 
is a persistant source, we expect that the disc is irradiated 
and therefore hotter than a standard steady state disc; note 
that an irradiated disc is also required to fit the UV contin- 
uum observed by Vrtilek et al. (1991). We therefore expect 
that the slope of the optical continuum will be considerably 
steeper (i.e. bluer) than that predicted for a steady state 
disc. Figure 2 illustrates that the observed optical spectrum 
is indeed much bluer than the simulated disc. (Note that the 
simulated spectrum has been arbitrarily normalized; thus 
this discussion concerns only the shape of the continuum, 
not the flux level.) However, it is unclear why the observed 
IR spectrum is shallower than the steady state disc model. 
If the disc is irradiated, we would expect the slope of the 
IR spectrum to be steeper than a standard disc. One pos- 
sibility is that the secondary star contributes to the flux 
of the IR continuum. A cool, late type (K/M) star would 
have a relatively shallow slope in the IR, whereas an earlier 
type (A/F) star would have a slope comparable to that ex- 
pected from a standard disc. In theory, therefore, flux from 
a cool secondary could cause the observed slope of the Sco 
X-l spectrum to be shallower than expected from the disc 
alone. However, if a late-type companion contributed suffi- 
cient flux to alter the slope of the continuum, that would 
imply a significant amount of the IR flux would be originat- 
ing from the secondary star. In this case, we would expect to 
see absorption features intrinsic to the secondary in the Sco 
X-l spectrum (see e.g. the modelling discussed in Paper I). 
As there is no evidence for such features in our spectrum, 
it seems unlikely that flux from the mass-donating star is 
causing the discrepancy between the observed and expected 
continuum slopes. It is therefore difficult to reconcile the 
continuum shape of the IR spectrum within the context of 
either a simple irradiated or steady state disc model. 

The differences between the Hel and Br 7 line strengths 
in our 1992 and 1997 spectra may indicate variability in the 
level of X-ray irradiation of the disc. After dereddening the 
1992 A'-band spectrum by Eb-v = 0.3, we find the slope of 
the continuum to be -3.34±0.08, i.e. somewhat steeper than 
the 1997 spectrum. The bluer continuum and the stronger 
emission lines in the 1992 spectrum are consistent with a 
higher level of X-ray irradiation than in the 1997 data. As 
we do not have any direct information about the X-ray state 
of Sco X-l during the 1992 observations, it is difficult to 
compare the X-ray behaviour of the source during the two 
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epochs quantitatively. However, such changes in the contin- 
uum slope are consistent with the behaviour expected from 
a disc which is being irradiated by a variable level of X- 
ray emission. Further, we note that the difference in the 
resolution of the two spectra and the highly changeable at- 
mospheric absorption in the region of the 2.059pim Hel line 
makes any conclusions derived from the apparently large 
change between the Hel line strengths at the two epochs 
uncertain. 

3.2 Sco X-2 

Our A'-band spectrum of Sco X-2 (GX349+2) was produced 
by combining our two observations for a total integration 
time of 80 minutes; the spectrum is presented in Figure 3. 
The prominent features are the Brackett 7 and Hel emis- 
sion lines from the accretion disc and/or heated face of the 
secondary. The presence of these lines directly confirm the 
identification of the IR counterpart to the X-ray source. The 
relatively low signal-to-noise (S/N ~ 8) despite the long ex- 
posure time is a result of the intrinsic faintness of the IR 
counterpart (K ~ 14.6). The line identifications, measured 
wavelengths, and equivalent widths are listed in Table 2. 

There is no evidence in our spectrum for the CO band- 
heads expected from a late-type secondary; however, the low 
S/N of our spectrum prevents us from drawing firm conclu- 
sions about the secondary's spectral type. It is interesting 
to note that Southwell, Casares, & Charles (1996; hereafter 
SCC96) found no evidence in optical spectra of Sco X-2 for 
the A6495 Cal/Fel feature, which is a signature of late G- 
and K-type stars (Home, Wade, & Szkody 1985). We also 
note our GX13+1 spectrum taken with CGS4 in 1995 did 
show clear evidence of the CO bandheads despite having 
a S/N comparable to that of our Sco X-2 spectrum shown 
here (Paper I). To date, the secondary in Sco X-2 has not 
been spectroscopically detected in either optical or IR wave- 
lengths. The lack of spectroscopic features associated with 
the secondary may indicate that the mass-donating star in 
this system is of an earlier type than the K/M secondaries 
detected in a number of LMXBs. Alternatively, the accre- 
tion disc contamination could be sufficiently high at both 
optical and IR wavelengths as to obscure late-type stellar 
spectral features in both our spectrum and the SCC96 op- 
tical spectrum. 

Several periodicities have been detected during observa- 
tions of Sco X-2. SCC96 suggested a period of ~14d, whereas 
Wachter (1997) found a period of 22.5±0.1h. As either of 
these periods could be the ~l-day alias of the other, it is 
unclear which of the two is the true orbital period (Barziv et 
al. 1997). The 14d period is clearly inconsistent with a main 
sequence secondary star; assuming a Roche-lobe filling sec- 
ondary, a spectral type of KOm would be required. However, 
in this case we again note that we might expect to see the 
prominent CO absorption features of a KOm secondary in 
our Sco X-2 spectrum. In contrast, a 22. 5h orbital period 
would be insufficient to contain a late-type giant, but also 
indicates that the companion star is evolved; in this case, 
the secondary would most likely be a subgiant. Several other 
XRBs with orbital periods ~20h have subgiant secondaries, 
including Sco X-l, Cen X-4 (P=15.6h; Shahbaz et al. 1993), 
and Aql X-l (P=19.1h; Shahbaz et al. 1996), suggesting that 
Sco X-2 may be similar. 



Table 5. Optimal subtraction of the GX13+1 companion star 



Spectral Type X 2 (DOF=146) / 



G5in 


495.0 


0.64±0.03 




ooy.o 




KOm 


421 3 


n A7-\-C\ fl9 


Klin 


647.6 


0.66±0.04 


K4m 


436.5 


0.46±0.02 


K5m 


284.5 


0.45±0.02 


MOm 


446.6 


0.52±0.02 


Mini 


424.7 


0.41±0.02 


M2m 


481.3 


0.48±0.02 


M4m 


465.8 


0.35±0.02 


G8iv 


714.2 


0.82±0.07 


KOiv 


507.6 


0.47±0.03 


K2iv 


529.2 


0.44±0.02 



3.3 GX13+1 

Figure 4 shows our K-band spectrum of GX13+1, produced 
by combining two spectra for a total on-source integration 
time of 88 min (S/N~ 30). The spectrum shows Br 7 and 
Hel emission, along with five 12 CO bands and three 13 CO 
bands which are characteristic of evolved late-type single 
stars (KH86). Marginal detections of Cal and Mgl are also 
present. The line identifications, measured wavelengths, and 
equivalent widths are listed in Table 2. 

3.3.1 The secondary star 

In Paper I we reported the detection of Br 7 emission in 
GX13+1, which confirmed the identity of the IR counter- 
part to the X-ray source. Also visible in the earlier spectrum 
were two CO bandheads; however, the much lower S/N of 
that spectrum only allowed us to constrain the secondary 
spectral type to be between K2 to M5. Together with an 
estimate of the distance, we also determined that the sec- 
ondary in this system must be a giant or subgiant. These 
constraints were consistent with previous estimates of the 
secondary's spectral type (e.g. Garcia et al. 1992). Using 
our new spectrum of GX13+1, we have now estimated the 
spectral type of the companion star in GX13+1 by the tech- 
nique of optimal subtraction, which minimizes the residuals 
after subtracting different template star spectra from the 
target spectrum. This method is sensitive to the fractional 
contribution of the companion star to the total flux /, where 
1-f is the "veiling factor" (Marsh, Robinson, & Wood 1994). 

First we determined the velocity shift of the spectrum of 
GX13+1 with respect to each template star spectrum by the 
method of cross-correlation (Tonry & Davis 1979). We then 
performed an optimal subtraction between each template 
star and the GX13+1 spectrum. The optimal subtraction 
routine minimizes the residual scatter between the target 
and template spectra by adjusting a constant /, which rep- 
resents the fraction of light contributed by the template star. 
The scatter is measured by carrying out the subtraction and 
then computing the x between the resultant spectrum and 
a smoothed version of itself. The constant / is therefore the 
fraction of light arising from the secondary star. The optimal 
values of / are obtained by minimizing \ 2 ■ 
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For this analysis we used a variety of giant and subgiant 
templates, ranging from G5 to M4. The optimal subtraction 
was performed in the spectral range 2.28-2.39 /jm in or- 
der to encompass the first four 12 CO bands, which are the 
most prominent absorption features in the GX13+1 spec- 
trum. The results of this analysis are presented in Table 5. 
The minimum \ 2 occurs for a K5iii companion star, where 
the secondary contributes about 45% of the flux at ~2.3 
/im. Figure 5 illustrates the result of the optimal subtrac- 
tion of the K5iii template from the GX13+1 spectrum. We 
note that the large IR variability which was observed by 
CN92 ruled out the possibility of the IR flux arising primar- 
ily from the mass-donating star; in other words, if the disc 
contribution to the IR flux was small, then we would expect 
the IR magnitude of GX13+1 to remain largely constant 
(as it does in GX1+4). The substantial IR variability (~ 1 
mag at K) indicates that we should expect X-ray heating 
of the disc and secondary to be a major contributor to the 
IR flux from the system. The result of our optimal subtrac- 
tion, which indicates that the K5iii secondary contributes 
less than half of the K-b&nd flux, is completely consistent 
with this expectation. 

The ratio of the equivalent widths of 12 CO to 13 CO 
depends upon luminosity class (see e.g. Dhillon & Marsh 
1995), ranging from 90 in main sequence stars to ~10 in 
giants (Campbell et al. 1990). To obtain a rough estimate 
of this ratio in GX13+1, we measured equivalent widths by 
first establishing the continuum as a linear function between 
marked points on either side of each feature. The flux was 
then determined by summing the pixels within the marked 
area and subtracting the continuum. Measurement of the 
ratio of the most clearly resolved 12 CO and 13 CO pair, the 
(2,0) bandheads, yielded a value of ~9, which is comparable 
to that expected for a field giant and inconsistent with a 
main sequence secondary. 

3.3.2 The P Cygni profile 

In addition to the absorption features from the secondary, 
the Br 7 line in our GX13+1 spectrum exhibits a prominent 
P Cygni profile. Similarly to Sco X-l, using the emission 
lines we have estimated the radial velocity of GX13+1 as 
well as measuring the location of the profile edges (see Table 
3). The inferred wind velocity, ~ 2400 km s -1 , is similar to 
that calculated for Sco X-l; as in Sco X-l, the outflow prob- 
ably originates in the disc rather than the mass-donating 
star. A high resolution spectrum of the Br 7 line is nec- 
essary to obtain accurate information about the outflow in 
this system. 

3.3.3 System parameters 

There have been several attempts at finding an orbital pe- 
riod in GX13+1. Groot et al. (1996) have observed a possible 
12.6±1 day modulation, while Corbet (1996) has reported a 
25. 2d period in the GX13+1 XTE ASM light curve. We note 
that the Groot et al. period is half of that found by Corbet. 
In contrast, Wachter (1996) found some evidence for a 19. 5h 
periodicity in Jf-band photometric data. Using Paczynski's 
formula for a Roche lobe filling star and an orbital period 
P = 25. 2d, the mean density for the secondary would be 



0.0003 gem -3 ; for a giant, this density corresponds to an 
approximate spectral type of K5. The K5iii spectral type 
we find for the secondary in GX13+1 therefore correlates 
well with the reported 25. 2d period, and supports the pos- 
sible identification of this period as an orbital modulation. 
In addition, the low Lx/L opt ratio in GX13+1, possibly the 
result of having a large region involved in X-ray reprocessing 
(i.e. a large disc), also reinforces the probability of a long or- 
bital period in this system (CN92). Our spectrum rules out 
an orbital period of 19. 5h, which would be insufficient to 
contain the orbital separation of a 1.4M Q neutron star and 
a Roche lobe filling K5m companion (M ~ 5Mq , R ~ 25R© ; 
Allen 1973). The nature of the 19. 5h periodicity is therefore 
unclear. 

Using an apparent K magnitude of 12 (CN92) and an 
absolute magnitude Mk = -3.8 for a K5iii star (Allen 1973; 
Koorneef 1983), together with a colour excess of Eb~v = 
5.7 (van Paradijs 1995), we find a distance to GX13+1 of 6.9 
kpc. There are several published estimates for the reddening 
(A v ) of GX13+1, ranging from ~17.6 to ~13.2 (see the dis- 
cussion in CN92), leading to an uncertainty in our distance 
calculation of approximately 1 kpc. Therefore we adopt a 
value of 7±1 kpc for the distance to GX13+1. This value is 
consistent with the previous estimate of 8-9 kpc, which was 
based upon the mean distance to the Galactic Centre and 
has been used generally as an estimate for the distances to 
all of the GBS (see e.g. Naylor & Podsiadlowski 1993). 

3.3.4 Evolutionary status 

On the basis of their X-ray properties, LMXBs have been 
divided into two classes, known as Z and atoll sources from 
the shape of their X-ray colour-colour diagrams (see van der 
Klis 1995 for a review). GX13+1 has been classified as an 
atoll (less luminous) source (Hasinger & van der Klis 1989, 
hereafter HK89), although earlier studies placed GX13+1 
into the category with the Z (high luminosity) sources (e.g. 
White, Stella, & Parmar 1988). HK89 suggested that the 
evolutionary history of the Z and atoll LMXBs could ac- 
count for their observational differences, with the Z sources 
having evolved secondaries and hence long orbital periods 
( > 15h) whereas the atoll sources would have main sequence 
secondaries and P or b ^ 5h. However, all of the proposed or- 
bital periods found to date in GX13+1 are substantially 
longer than the ~5h periods predicted for an atoll source 
and in fact are similar to the periods found in the Z sources 
Sco X-l and Cyg X-2. Additionally, our spectrum shows un- 
equivocally that the mass-donating companion in GX13+1 
is an evolved late-type star, making a ~5h orbital period 
impossible. 

There are several possibilities for resolving the discrep- 
ancy between the observed nature of GX13+1 and the pre- 
dictions arising from the Z/atoll classification scheme. First, 
GX13+1 may be mis-classified as an atoll source. Both the 
high X-ray luminosity of GX13+1 and the detection of ra- 
dio flux from the source are more typical of the luminous 
Z LMXBs than atoll sources; additionally, a ~60Hz QPO 
has recently been found (Homan et al. 1998). However, its 
colour-colour diagram is not consistent with the X-ray prop- 
erties of the six known Z sources. A second possibility is 
that the existence of "hybrid" systems such as GX13+1, 
GX9+9, and GX9+1, which exhibit both Z and atoll char- 
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acteristics, should be considered as a third "intermediate" 
class of objects where the Z and atoll classes overlap (e.g. 
Kuulkers & van der Klis 1998). Finally, it is possible that the 
evolutionary distinction between Z and atoll LMXBs does 
not arise from a straightforward dichotomy between evolved 
and main sequence secondaries. For example, some systems 
could have evolved companions but lower accretion rates 
and/or weaker magnetic fields than canonical Z sources. In 
this case, sources could exhibit atoll-type X-ray properties 
as well as long periods; in addition to GX13+1, the atoll 
source AC 211 (P or b = 17. lh) could fit into such a category 
(van der Klis 1992). It would then become difficult to make 
predictions about the orbital period and nature of the com- 
panion star in a given system purely on the basis of its X-ray 
behaviour. However, we note that the Z/atoll evolutionary 
dichotomy could still arise from fundamental differences in 
the types of mass-donating star in the two classes, but in a 
more specific manner than a simple division between evolved 
and main sequence companions. 

3.4 GX5-1, GX17+2, 4U2129+47 

In Figures 6 and 7 we show the spectra obtained of candidate 
IR counterparts to the GBS GX5-1 and GX17+2 and the 
LMXB 4U2129+47. To date, no IR counterpart has been 
confirmed for any of these elusive LMXBs. 

34.1 GX5-1 

Other than Sco X-l, GX5-1 is the brightest of the persistent 
LMXBs, but no X-ray periodicity has been found nor has an 
IR counterpart been identified (van der Klis et al. 1991). An 
Einstein X-ray position and a precise radio position exist for 
GX5-1, but candidates for the IR counterpart to this heav- 
ily obscured source only become visible at K (see e.g. the 
finding charts in NCL91). We obtained spectra of the two 
stars closest to the radio error circle, designated by NCL91 
as stars 502 and 503. Neither of the spectra, shown in Fig- 
ure 6, show any evidence for the signature Br 7 emission we 
would expect in an LMXB. Both stars 502 and 503 are faint, 
with K magnitudes of ~14 and ~15 respectively (NCL91); 
combined with the relatively short integration time (16 min 
on each source) , the generally poor quality of the two spec- 
tra (S/N ~ 3) is unsurprising. However, in such a bright 
X-ray source, we might expect the emission lines from the 
disc to be especially strong; while not conclusive, our spectra 
therefore cast doubt on the potential of either of these two 
stars to be the IR counterpart of the X-ray source. We note 
that there is a star northeast of star 502 (designated as star 
513 by NCL91) which also lies within the radio error circle 
of GX5-1 and as such is also a prime candidate for the IR 
counterpart. Due to the star's extreme faintness (K ^16), 
we did not obtain a spectrum for this candidate. Obtain- 
ing photometric and spectroscopic information on star 513 
should therefore be a primary objective for future IR obser- 
vations of GX5-1. 

3.4.2 GX17+2 

On the basis of its X-ray position, the GBS GX17+2 was 
optically identified with a G star known as "star TR" more 



than 25 years ago (Tarenghi & Reina 1972; Davidsen, Ma- 
lina, & Bowyer 1976). The location of GX17+2 was fur- 
ther refined by a subarcsecond radio position (Hjellming 
1978) but the putative counterpart has shown almost no 
optical variability (e.g. Margon 1978). We note, however, 
that Deutsch et al. (1996) have found a small discrepancy 
between the optical and radio positions of GX17+2. The 
IR counterpart to star TR (K ~14) has shown variability, 
but a consistent fit for the colours, extinction, distance, and 
spectral type could not be found (NCL91). Based on the 
variety of X-ray periods which have been suggested for this 
source, Bailyn & Grindlay (1987) theorized that GX17+2 is 
a triple system, with a giant star in orbit around a short- 
period LMXB. However, the discrepancy between the ex- 
tinctions derived from optical and X-ray measurements led 
NCL91 to conclude that star TR is a foreground star unas- 
sociated with the X-ray source which is superimposed upon 
a highly reddened IR counterpart to GX17+2. Our K-baxid 
spectrum, shown in the top panel of Figure 7, shows no ev- 
idence either for emission lines or absorption features from 
a late-type IR counterpart obscured by star TR. Indeed, 
the featureless spectrum is consistent with an early G-type 
evolved star or a main sequence star earlier than G5 (KH86), 
as would be expected from star TR if it is an unassociated, 
foreground G star. Therefore, our spectrum has too short 
an integration time and/or insufficient resolution to show 
features from a faint, variable IR counterpart to GX17+2 
which may be partially hidden behind star TR. 

3.4.3 4U2129+47 

The neutron star LMXB 4U2129+47 is known to have an 
orbital period of 5.2h and has an optical counterpart, V1727 
Cyg (Thorstensen et al. 1979). Radial velocity studies pro- 
duced estimates for the masses of the compact object and 
secondary star of O.6±O.2M0 and 0.4±0.2Mq respectively 
(Thorstensen & Charles 1982). Then, after entering X-ray 
quiescence in 1983, all evidence of photometric variability 
and radial velocity variations disappeared (e.g. Garcia et al. 
1989, hereafter G89). The optical counterpart visible dur- 
ing the low state is an F8iv (Cowley & Schmidtke 1990) 
whose colours were initially observed to be inconsistent with 
a normal star (Kaluzny 1988). However, subsequent obser- 
vations failed to confirm this abnormality in the optical 
colours of the F star (Deutsch et al. 1996). G89 proposed 
that 4U2129+47 is a triple system, although it has also been 
suggested that the F star is a foreground source unassoci- 
ated with the system. Our _K"-band spectrum of the F star 
appears in the bottom panel of Figure 7. The spectrum is 
featureless (the spike at ~2.05 /jm is a residual of the tel- 
luric absorption feature removal); we see no indication of 
emission from an accretion disc, nor do we see any evidence 
for absorption features from the Mv star suggested by G89 
as a hypothetical third member of the triple system. How- 
ever, we note that even in the K-han&, the flux from a faint 
Mv star would likely be obscured by the brighter F star. 



4 CONCLUSIONS 

We have presented JT-band spectra of the LMXBs Sco X-l, 
Sco X-2, and GX13+1. The IR spectrum of Sco X-l exhibits 
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the emission features expected from a luminous accretion 
disc but does not show any absorption features arising from 
the system secondary, reinforcing our conclusion that the 
secondary is a subgiant of spectral type G5 or earlier. The 
spectrum of the proposed counterpart to Sco X-2 shows Br 
7 emission from the accretion disc, confirming the identifica- 
tion. The if-band spectrum of GX13+1 exhibits both Br 7 
emission and CO absorption bands; optimal subtraction in- 
dicates that the secondary in GX13+1 is most likely a K5ni, 
with the accretion disc contributing ~50% of the K-band 
flux. In addition, the Br 7 emission lines in both Sco X-l 
and GX13+1 show P Cygni profiles, indicating the presence 
of outflows in both systems with velocities ~2000kms _1 . 
Spectra of two stars in the error circle of GX5-1 do not 
show the Br 7 emission signature which would identify the 
IR counterpart. Similarly, spectra of the G star at the posi- 
tion of GX17+2 and the F star at the position of 4U2129+47 
fail to exhibit any evidence that they are physically associ- 
ated with the X-ray sources. 

It is interesting to note that while our GX13+1 spec- 
trum clearly shows the features of the K giant secondary, 
neither Sco X-l nor Sco X-2 appear to exhibit any spectral 
lines characteristic of late-type (late G to M) stars. Signif- 
icant similarities in the behaviour of Sco X-l and Sco X-2 
have also been seen in the optical and X-rays, indicating 
a fundamental similarity in the nature of the two systems. 
Furthermore, the clear differences in these spectra occur de- 
spite the fact that GX13+1, Sco X-l, and Sco X-2 all have 
relatively long orbital periods and are luminous, persistent 
XRBs. Far from supporting the idea that most of the GBS 
have late-type evolved secondaries, this result leads us to 
speculate that the mass-donating stars in Sco X-l and Sco 
X-2 may prove to be very different than previously assumed. 

There have been two primary schemes proposed for the 
classification of LMXBs. In a statistical study of a flux- 
limited sample of LMXBs, Naylor and Podsiadlowski (1993, 
hereafter NP93) asserted that the GBS distribution is as- 
sociated with that of M giants in the Galactic bulge. They 
classified Sco X-l, Sco X-2, and GX13+1 (along with seven 
other XRBs) as bulge sources on the basis of their location, 
the characteristics of their X-ray spectra, and the fact that 
they are persistent X-ray sources. In contrast, the Z/atoll 
scheme categorizes LMXBs based on the shape of their X-ray 
colour-colour diagrams (HK89). With this method, six of the 
sources classified as GBS by NP93, including Sco X-l and 
Sco X-2, are labelled as ^-sources, whereas the remaining 
GBS, including GX13+1, are placed in the atoll category. 
Note that the Z and atoll sources appear to be randomly 
scattered amongst the bulge and disc populations catego- 
rized by NP93. HK89 suggested that a primary difference be- 
tween the Z and atoll sources may originate with the mass- 
donating stars; namely, that the Z sources have evolved sec- 
ondaries while the atoll sources have main sequence compan- 
ions. It is clear from our spectra that the dichotomy between 
Z and atoll sources cannot be this straightforward, as both 
Sco X-l and GX13+1 have evolved secondaries. However, 
if the spectral type of the secondaries in Sco X-l and Sco 
X-2 prove to be considerably earlier than the K giant found 
in GX13+1, then the evolutionary scenarios of Sco X-l and 
Sco X-2 are indeed likely to be quite distinct from sources 
such as GX13+1. We also note that the Z source Cyg X-2 
has an evolved A9 companion (Casares, Charles & Kuulkers 



1998); Sco X-l also has an evolved secondary of type earlier 
than G5, and Sco X-2 may be similar. 

We therefore hypothesize that the secondaries in the six 
Z sources may all be early-type evolved stars. We speculate 
further that the distinction between the evolutionary his- 
tory of the Z and atoll sources may not be that the former 
have evolved secondaries and the latter have main sequence 
companions, but instead that Z sources have evolved early- 
type companions while most other LMXBs have late type 
(either evolved or main sequence) secondaries. If this were 
the case, it may also help to explain why the IR counter- 
parts for the other three Z sources (GX5-1, GX17+2, and 
GX340+0) have proven so elusive (although GX17+2 is a 
special case due to the foreground star). Until now it has 
been assumed that the secondaries in these systems would 
be late-type evolved stars and therefore bright in the K band 
(V — K ~ 4). Yet the few candidate stars found in deep K- 
band images of the GX5-1 and GX340+0 fields have failed to 
show photometric or spectroscopic evidence that they are in 
fact the counterparts. To elude detection, a late-type evolved 
counterpart would have to be at a greater distance than we 
believe from estimates of the X-ray column densities. How- 
ever, if the secondaries in these systems are early-type stars 
(with V — K ~ 0), then the counterparts would naturally 
be fairly faint in the K-h&n&. It would then be unnecessary 
to invoke a significantly larger extinction and/or distance to 
these sources in order to explain our inability to locate their 
IR counterparts. 

Our spectra provide clear evidence that neither the 
bulge/disc nor the Z/atoll schemes can be explained by a 
simple distinction between evolved and main sequence sec- 
ondary stars. However, we note that the early- type/late- type 
dichotomy proposed here could fit into the Z/atoll classifi- 
cation scheme without having to disregard the bulge/disc 
association made by NP93. Although currently the evidence 
for this claim is limited, we believe that this hypothesis is 
well worth exploring. 
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Figure 1. _ft"-band spectrum of Sco X-l (K ~ 11.8). Strong emission lines of He I, He II, and Br 7 are evident, but no absorption 
features from the secondary can be distinguished. Note that the Br 7 line has a P Cygni profile. 
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Figure 2. The continuum of the observed optical (thick line) and ii"-band spectra of Sco X-l are plotted together with the calculated 
spectrum (F\) of a steady state accretion disc with T ou t = 4200 K (dashed line) for comparison. The parameters used for calculating 
the temperature were an outer disc radius of 6.3 X 10 10 cm (Vrtilek et aJ. 1991), L x = 3 x 10 37 ergs-^Bradt & McClintock 1983), and 
the inner radius of the disc at the last stable orbit around the neutron star. The simulated disc spectrum has been arbitrarily normalized 
to match the flux of the optical spectrum at 0.9/im. 
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Figure 3. iiT-band spectrum of the GBS Sco X-2 (K ~ 14.6). Br 7 and He I emission lines are clearly visible, confirming the identity 
of the IR counterpart to the X-ray source. 
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Figure 4. A"-band spectrum of the GBS GX13+1 (K ~ 12.0). Visible features include Br 7 and He I in emission and absorption 
bands of 12 CO and 13 CO, and with marginal detections of Cal and Mgl. The Br 7 line has a P Cygni profile with a strong absorption 
component. 
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Figure 5. This figure shows the spectrum of GX13+1 and that of the template K5m star, as well as the residual after optimal 
subtraction. The region including the first four 12 CO bandheads was used to perform the subtraction. Also shown is the A- type star 
used for removal of the telluric absorption features. The flux level has been normalized and shifted for ease of plotting. 
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Figure 6. The spectra of two candidates for the counterpart to GX5-1, stars 502 (top) and 503 (NCL91). 



© 0000 RAS, MNRAS 000, 000-000 



12 R.M. Bandyopadhyay et al. 




2 2.1 2.2 2.3 2.4 

wavelength (jum) 



Figure 7. The top panel is the spectrum of "star TR", the G star at the position of GX17+2. The lower panel is the spectrum of the 
F star at the posision of 4U2129+47. (The "feature" at 2.05 (im is a residual of the telluric absorption line removal.) 
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